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Universality in s-wave and higher partial wave Feshbach resonances: an illustration 

with a single atom near two scattering centers 

Shangguo ZhiQ and Shina Tai|i] 
School of Physics, Georgia Institute of Technology, Atlanta, Georgia 30332, USA 

(Dated: May 31, 2013) 

It is well-known that cold atoms near s-wave Feshbach resonances have universal properties that 
are insensitive to the short-range details of the interaction. What is less known is that atoms near 
higher partial wave Feshbach resonances also have remarkable universal properties. We illustrate 
this with a single atom interacting resonantly with two fixed static centers. At a Feshbach resonance 
point with orbital angular momentum L > 1, we find 2L -\- 1 shallow bound states whose energies 
behave like 1/R '^^ when the distance R between the two centers is large. We then compute 
corrections to the binding energies due to other parameters in the effective range expansions. For 
completeness we also compute the binding energies near s-wave Feshbach resonances, taking into 
account the corrections. Afterwards we turn to the bound states at large but finite scattering 
volumes. For p-wave and higher partial wave resonances, we derive a simple formula for the energies 
in terms of a parameter called "proximity parameter" . These results are applicable to a free atom 
interacting resonantly with two atoms that are localized to two lattice sites of an optical lattice, and 
to one light atom interacting with two heavy ones in free space. Modifications of the low energy 
physics due to the long range Van der Waals potential are also discussed. 

PACS numbers: 67.85.-d,31.15.ac,03.65.Ge 
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I. INTRODUCTION 

Cold atoms with de Broglie wave lengths that greatly 
exceed the range of inter- atomic potential are known to 
exhibit universal properties that are mainly determined 
by the s-wave scattering length hq. Interactions in higher 
partial wave channels are usually suppressed by the cen- 
trifugal barrier, unless the atoms are near a p-wave [iHlJl , 
d-wave [14| - [lq . or higher partial wave jl7H2(j | Feshbach 
resonances. 

The universaUty is particularly remarkable when the 
atoms are near a broad [2l|, [22| s-wave Feshbach reso- 
nance, near which oq is large but the effective range tq 
is small [2j, [2J|. For example, if oq is positive and large 
there is a shallow two-body bound state with energy en- 
tirely determined by an rather than the full atomic details 
of the interaction [2^, [2^] : 



E^. 



2fiaQ ' 



where /x is the reduced mass, and h is the Planck con- 
stant over 27r. For three or more atoms with large scat- 
tering length, there is another remarkable manifestation 
of universality: the Efimov effect [23, [23; i'^ particular, 
at oo = ±oo, there are an infinite sequence of three-body 
shallow bound states whose energies form a geometric 
sequence with a ratio that is determined by the quan- 
tum statistics and the mass ratios of atoms j26, 29, 30], 
but not affected by the details of short-range interac- 
tion |26l - [28| . There have been predictions of analogous 



effects in two-atom systems with large s-wave scattering 
length [H. 

Unlike the s-wave resonance, higher partial wave Fesh- 
bach resonances have not been extensively analyzed from 
the perspective of universality. There have been efforts 
to look for Efimov effect near p-wave resonances |32h34| | , 
but Nishida pointed out that Efimov effect is not possible 
in such a case |35| . 

In this work, we illustrate that near higher partial wave 
resonances, although there is no Efimov effect in three di- 
mensions, the concept of universality remains powerful. 
By universality we mean that there are many physical 
properties controlled by one or two effective parameters 
for the interaction, rather than by the full atomic details. 
Consequently different atomic species near higher partial 
wave Feshbach resonances can exhibit very similar be- 
haviors. 

We study a single free atom of mass M having short- 
range interactions with two identical static centers. The 
interaction is tuned to a resonance in the Lth partial wave 
channel. Here L is the quantum number for the relative 
orbital angular momentum between the atom and one 
center. If the L-wavc scattering volume a^ is infinity, 
we find that there are 2L + 1 shallow bound states with 
energies entirely controlled by the L-wave effective range 



E= - 



XLr 



2M {-rL)R^^+^' 



m = —L, —L + 1, 



,,L (1) 



if i > 1. Here R is the distance between the two centers, 
m is the projection of orbital angular momentum along 
the line connecting the two centers. 



* Electronic address; 
T Electronic address: 



shangguo.zhu@gatech.edu] 



shina.tan®physics.gatech.edu| 



Xi„^2(2L + l)!!(2L-l)!! 



2L 



(2) 



and 



2L 



(2L)\ 



arc binomial coefficients. 



L + m J (L+m)!(L-rn)! 

aL and r^ are defined as parameters in the low-energy 
effective range expansion for the scattering phase shifts, 
Eq. ((131). Note that tl is negative ff L > 1 j4, 36_ sl]. 
We shall also study the bound states when a^ is large 
but finite in Sees. El El and ED 

Even when a^ = ±00, there are corrections to Eq. ([1]) 
due to a nonzero radius of the interaction. Such correc- 
tions decay more rapidly at large R than the dominant 
term, so that at large values of R Eq. ([T|) becomes increas- 
ingly more accurate. We determine these corrections in 
Secinil 

At the s-wave resonance (ao = ±00), the problem of a 
single atom interacting with two centers has been studied 
in Refs. [H HO, IMSl. It was found that the shallow 
binding energy 



E^- 



'Xo 



2MR^ 



(3) 



and is independent of the effective range tq at the lead- 
ing order ^ [s^, ISBl- Here xo ~ 0.567143 is the 
real solution to the equation xo = e^^°- This inverse 
square dependence of the energy on R qualitatively ex- 
plains the presence of Efimov effect for two heavy atoms 
interacting resonantly with a light atom. Within the 
Born-Oppenheimer approximation the two heavy atoms 
arc effectively interacting via a potential equal to the 
binding energy for the light atom. Since the potential 
decays like 1/R^ and is strong compared to the kinetic 
energy of the heavy atoms, there are a sequence of three- 
atom bound states with energies forming a geometric se- 
quence [2i,[33,li3- 

For two heavy atoms interacting with a light atom near 
higher partial wave Feshbach resonances, Eq. ([T]) shows 
that Efimov effect is absent, in agreement with Ref. |35| . 
since the effective potential induced by the light atom 
now decays like l/R^ for p-wave resonance, 1/R^ for d- 
wave resonance, etc. 

Despite the absence of Efimov effect, the simple power 
laws for the shallow bound states in Eq. ([!]) illustrates the 
elegance of universality in higher partial wave resonances. 

In Sec. lVII we derive a more general simple approximate 
formula, Eq. ([35]), when laLl^/'^^^^^' and R are both large 
and finite. 

One could experimentally confirm Eq. ([Ij and Eq. (|38p 
using a light atom and two heavy atoms, or using a 
free atom having confinement induced scattering reso- 
nance [43, [4^ with two other atoms pinned to two lattice 
sites of an optical lattice. 

The outline of the rest of this paper is as follows. 

In Sec. |TT]we derive the general equations for the shal- 
low bound states of an atom near two centers. 

In Sec. mil we derive large- i? expansions for the bound 
states at resonance (a^ = 00). 

In Sec. HV] we derive the 0{l/aL) corrections to the 
wave numbers of the bound states slightly off resonance. 



In Sec. El wc compute the binding energies approxi- 
mately when |aL|^/'-^^+^' and R are both large and may 
be comparable. 

In Sec. IVII we derive a simple formula for the binding 
energies near a p-wave or higher partial wave resonance. 

In Sec. lVIII we consider the effects of a long range Van 
der Waals interaction between the atom and each center. 

In Sec. IVIIII we summarize our findings, discuss some 
technical limitations of this work, and argue that despite 
such limitations, our results will be useful in various ex- 
perimental scenarios. 



II. GENERAL PROPERTIES OF THE BOUND 

STATES 

We consider an atom of mass M attracted by two fixed 
identical scattering centers that are separated by distance 
R and located along the z axis, with Cartesian coordi- 
nates si = (0,0,-i?/2) and 83 = {0,0,+R/2). The po- 
tential experienced by the atom is assumed to vanish un- 
less the atom is very close to one center. This idealized 
model is a good approximation for neutral atoms, which 
have a long range Van der Waals potential that decays 
rapidly at large distance; effects of this long range poten- 
tial are discussed in Sec. IVIII 

Each scattering center is characterized by phase shifts 
di{k) where I = 0,1,2, . . . is the orbital angular momen- 
tum of the atom about the center and k = in is the wave 
number. The energy of the atom is 



E^ 



2M 



2M 



(4) 



For a bound state we take k > 0. 

Due to the axial symmetry of the system about the z 
axis and the parity symmetry about the origin, we may 
consider bound states with a certain orbital magnetic 
quantum number m and a certain parity a £ {-|-1,— 1}. 
Outside the range of potential the wave function of the 
atom is 

00 

+ a{-iyh['\tKr2)Yi^ir2)\, (5) 



which satisfies 



■0mcr(-r) = aipmair). 



(6) 



Here r is the position vector of the atom, r^ = r — Sq 
is the position vector of the atom relative to a center, 

fc = Va/ra, Yim is thc sphcrical harmonic, and h^ (x) 
is the spherical Hankel function of the first kind. The 
coefficients C™°' and the parameter k should be chosen 
such that in the vicinity of a center located at s^ (a = 



1, 2) the wave function takes the form [25| (O at one center only, such as the one located at Si. At 

ri < R we have a useful expansion: 

oo 

^™<,(r)- ^ Qr[ji(fc^a)cotJ;(fc)-n,(fcr„)]yi„(f„), 

l=\m\ °° 

(7) h'i'\iKr2)YUr2) = {-iy ^ Fl™'M)jH*'^^i)>l'™(fi), 
where j; (x) and n; (x) are the spherical Bessel functions of i'=\m\ 

the first kind and the second kind, respectively. Because (8) 

of the well-defined parity, we need to impose condition where 



FiTix) 



inin{l,l )- 

E 

3=0 



j\{2j + 2m- 1)!! f I 
(.? + 2m)! 




r 



(1) 



l)i2l' + !)(/ + m)\il' + m)\ h\^i,_^_^{ix) 



{I - m)\{l' - m)\ 



i-ix) 



m+j 



(9) 



and 



f-(n~jV- 



is the binomial coefficient. Substi- 



.(1 



tuting the identities h^ (x) = ji {x) + irii {x) and (|5]) into 
Eq. ©, and noting that Fllf{x) = F™(.t), wc find 



CXD 



OO 



^™.(r)= >J {[cr+ >: 1CrFP{nR)\M^Kn) 



l=\7n\ 



l' = \m\ 



+ iCrniiiKn)}Yi,rr{h). (10) 



Comparing the above formula with the required form of 
the expansion of the wave function in Eq. ([7]), we find a 
set of linear equations for the unknown coefficients C™"': 



J2 K]^.^(n,R)C]? 

V = \rn\ 



0, 



l>\ml (If) 



where 

KI^,''{k,R) = [l + icotSi{iK)]Sw + (tFIP^{kR). (f2) 

For each pair of values {m,a), we must adjust k in or- 
der to get a nonvanishing set of coefficients C""^. Each 
positive value of k corresponds to a bound state, whose 
energy is given by Eq. Q. For shallow bound states 
with k ~ IK close to zero, we shall use the effective range 
expansion for the phase shifts 5i [4^, |4^ : 



U2/+1 



cot5,(fc) = -l+r,^+r;^+rr^+0(fc«). (13) 



The coefficients r[, r" etc arc known as shape parameters. 
We will use Eqs. ([TT]) ([12]) ^^ to derive well-controlled 
expansions for many shallow bound states in Sec. IIIII 

At a given m, the two parities a = ±1 give rise to two 
different states. We can tell which one has a lower energy 
by a simple symmetry analysis. It is easy to see that 



i>niaix,y,-z) = Cr(-l)"Vmcr(a;, y,z). 



(14) 



If tT(— 1)™ = —1, the wave function vanishes on the xy 
plane, and the energy is the same as we would obtain 
by imposing a hard wall potential there; thus the closer 
the two centers, the closer the hard wall is to one center, 
and the higher the energy, hence the energy increases 
as we reduce the distance between the two centers. In 
the language of the chemical bond theory, the state with 
parity a ~ (—1)™+^ may be called an anti-bonding state. 
We have also verified numerically that for a = (—1)™, 
for which the wave function is an even function of z, the 
energy decreases as we reduce the distance between the 
two centers; see Sec. E) Thus 



(-!)'", 
(-f)'"+\ 



for bonding states, 
for anti-bonding states. 



(15) 



For a particular m, depending on the interaction between 
the atom and the centers and the distance R between the 
two centers, the atom may have two shallow bound states 
one from each parity, or a single shallow bound state with 
CT = (— 1)™, or no shallow bound state. More details will 
be shown in Sec. Ivl 



III. BOUND STATES AT RESONANCE 

In this section we study the shallow bound states at 
an L-wave scattering resonance. By this resonance we 
mean that a^ diverges but the scattering volumes ai for 
all other partial waves (with / 7^ L) remain finite. We 
shall call the partial wave channel with orbital angular 
momentum L the resonant channel, and all other chan- 
nels non-resonant channels. 

At large R we find 2i -t- 1 shallow bound states with or- 
bital magnetic quantum numbers m = — L, —L + 1, . . . , L 
and parities a = (-1)™. Using Eqs. (HH dH]) ^, we 
can expand k{R) for these states in powers of l/R. It 
is important to bear in mind that Eq. (jf f p contains the 
effects of all partial wave channels, not just the resonant 
channel. 



At 



an s-wavG resonance, 



0.567143 



R 



0.180948ro 0.0688446r2 -2.66638ai + 0.0332811r3 - 0.00485019r^ 



R 



i?2 



iJ3 



-0.865928airo + COlBUGTr* - 0.00544592ror;, 

(-45.655602 - 0.473306air^ + 0.0106396r^ - 0.00483849r§r[, + 0.0000520023r^' + 0.428823a?ri) 

R? 



O 



Vi?6 



(16) 



r 



The leading order term in the above series is consis- 
tent with previous theoretical calculations |26l . [30ll40l - t42| . 
The binding energy E has the — 1/i?^ dependence, which 
leads to the emergence of three-body Efimov states. 
Moreover, the effects of the parameters rg, Tq, Tq etc 
in the effective range expansion Eq. (J13p appear only in 



higher order terms. Also, the effects of other partial wave 
channels are all higher order corrections. In this sense, 
the resonant partial wave channel is dominant, in agree- 
ment with our intuition. 

The results at p, d, /-wave resonances are listed in the 
following. 



J 



1. For p-wave resonance: 
(a) m = 0, 

^ r. an 288 - 96aori -f Ta^jf ^3 (36 - 20aori 4- 5agr?) 



~i?3/2 



4i? 



+ 



32r2i?2 5(-ri)5/2i?5/2 

1536 - 864aori + 288agr^ - 25agr? - Mrir{ \/3 (180 + 924aori - 595agr2 + TOaj^rf ) 



O 



(i.) 



128rfi?3 



35(-ri)V2i?7/2 



(17) 



(b) m = ±1, 



a/6 



1- 



2^/6 



9 



57\/6 



1083 4- 4riri 



' \/^7Ti?3/2 

2. For d-wave resonance 
(a) TO = 0, 



2rii? (-ri)3/2i?3/2 8^2^2 5(_^^)5/2^.5/2 16^3^3 



1521^/6 -nf— "l 

"70(-ri)7/2i?7/2+ (,^4 J 



(18) 



/540 



^7ii?5/2 



1- 



Oq 



23a§ 



12i? 288i?2 



265ao _ 9aa i 15 
3456 4 "*" ra 

i?3 



15a^ 



-7-2 



i?V2 



165888 



16 



65ao 
4r2 



i?-4 



o 



i?9/2 



(19) 



(b) TO = ±1, 



/360 



-?'2 



i?5/2 



^ + 77 15r^ 27 (1600 - 720air2 + 13afri) 3937503 



i?3 



128r^i?f5 



45 (36afrir| -f 2G0r^ - 9air2r0 / 1 

8^p8 +^1;r9 



4i?7 



(20) 



(c) TO = ±2, 



90 



-^2 



i?5/2 



15 



15r(, 



2025 7875a3 



540V10 



1125r; 



2r2-R3 4r2i?5 8^2^6 2i?7 (_j.2)5/2^15/2 8^3^8 



O 



i?9 



(21) 



3. For /-wave resonance: 
(a) m = 0, 



1- 




flo 



40i? 3200i?2 



7Q„2 3121a;j 6a]_ 19731a;S STapai q,/7n „2 

'■^"O 128000 "^ 5 I 819200 ^ 100 OV'UUo 



i?3 



i?4 



;il'9/2 



3ai 15(147 -lOasrg) Slaj TSTSrJ, 9ai(-1407 + SOaarg) 



5i?3 



2r3-R5 



50i?6 4r2i?7 



2r3i?8 



O 



(i) 



o 



ii) 



(22) 



(23) 



V18900 



15(168-502^3) 1575r!, 4862025a4 225 (155232 - 11760a2r3 + 95a; 



36673560005 



R 



11 



O 



2r2i?7 



R 



12 



4^9 



32r|i?io 



2 2\ 
2'3J 



(24) 



(d) TO = ±3, 



V3T50 



315 525r!, 694575a4 628425 137525850a5 



2r3i?5 4r2i?7 



2i?9 



8r2i?,io 



i?i 



O 



i?i 



(25) 
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From Eq. P7|) - (P5|) . we see that at the L-wave reso- 
nance, the leading order term in the expansion of k is 
proportional to i}(2L+i)/.! - The first order correc- 

tion is due to the partial wave channels with orbital an- 
gular momentum quantum numbers equal to L or \m\, 
while the higher order corrections arc also influenced by 
other partial wave channels. By expanding k to high 
powers of 1/R, one will see the effects of all the effec- 
tive range expansion parameters in all the channels that 
satisfy I > \m\. 

In general, the leading order terms for the binding en- 
ergies at the L-wave resonance are 



E = 



hf_Xo 

2M {-rL)R'^^ + ^ 



L = 0, 
L> 1, 



(26) 



where xo ~ 0.567143 is the real solution to the equation 
Xo = e^^" , —L < m < L, and 



XLt. 



2(2L+1)!!(2L-1)!!(2L)! 

(L- to)!(L + to)! 



(27) 



To obtain Eq. (P5|) at L > 1, we used the following ex- 
pansion at small x: 



Fl\{x)^{-1) 



L— m+l ALr; 



2x 



2L + 1 



0{x 



-2L\ 



(28) 



When two heavy particles interact resonantly with a 
light particle by short-range potential, within the Born- 
Oppenhcimcr approximation they experience an effective 
potential equal to the binding energy E for the light par- 
ticle. At an s-wave resonance, the —\/R^ dependence of 



the effective potential results in the appearance of three- 
body Efimov states [40|. At higher partial wave reso- 
nances [L > 1), however, the three-body Efimov effect 
seems unlikely, because the Schrodingcr equation for two 
heavy particles with effective potential ex — l/i?2^+i has 
no scaling symmetry. 



IV. CORRECTIONS TO k SLIGHTLY OFF 
RESONANCE 

Since experimentally the scattering volume a^ can not 
be tuned exactly to infinity, we would like to know the 
corrections to the expansions in Eqs. p^ -(P5 |) when ql 
is large but finite. We use Eqs. ^ ^ i^^ to calculate 
these corrections to the leading order in l/a^. Our re- 
sults are listed in Table HI In general, in the very vicinity 
of a resonance 



I resonance 



K + , 

I resonance \/—X 



(l+Xo)ao ' 

^L + l/2 



L = 0, 
L> 1. 



(29) 



To guarantee the applicability of the expansions in 
Eqs. (|T6|) - (|25l) . the above corrections must be small com- 
pared to the leading order term for k. This condition is 
met if |aL|^/'-2-^+^) ^ R. Moreover, the effective range 
expansion Eq. ([T^ requires \kre\ <C 1, where r^ is the ra- 
dius of the short-range potential between the atom and a 
scattering center. In many systems rg can be character- 
ized by the effective range r^, namely Tg ^ jrij^/'^'^i+i)^ 
For such systems we obtain the domain of applicability 



m = 
m = ±l 
m = ±2 
m = ±3 



R 



3/2 



ii 



5/2 



R 



7/2 



2\/3\/ — ^l^l 6\/15-^ — r2 a2 30^ 70^ — r^^a^ 

fl3/2 flV^ fl^/^ 



Vo-^/— r^aj^ 6\/10y/ — r2O,2 15\/2T0^/^^r^a^ 

^5/2 ^7/2 



3\/K)V^rJa2 30V2Ts/^r^a3 

rV^ 



TABLE I: The 0(l/a_t) corrections to k slightly away from 
resonance. 



of the expansions in Eqs. (|T6 |) - (f25l) : 

|rL|i/(-2^+i)«i?«|aL|i/(2L+i). (30) 



NUMERICAL RESULTS FOR THE BINDING 
ENERGIES NEAR A RESONANCE 



resonance, Eq. ([32|) shows that the bound states with 
different values of |to| disappear at different critical dis- 
tances. The binding energy vanishes either quadratically 
or linearly when R approaches Re from below, depending 
on L: 



E: 



(|ao|-i?)^ 

Mfl,|a,,r^|(-Rc-i?), 



(2L + l)h 



L = 0, 
L> 1. 



(33) 



If cl is large and positive, we have either one or two 
shallow bound states for each orbital magnetic quantum 
number m ranging from —L to i, depending on the dis- 
tance R between the two centers. For all R there is a 
"bonding state" with parity cr = (— l)™. liR> Re where 
Re is given by Eq. ([5^ there is an extra shallow state with 
parity a = (—1)"'+^, namely the "anti-bonding" state 
discussed in Sec. |lll When the two scattering centers are 
far apart [R — ^ oo), the wave functions of both shallow 
states are localized near two separate centers, and their 
energies both exponentially approach the binding energy 
due to a single center: 



When R and jaij^/'^'^'''^'' are both much larger than 
the radius r^ of the potential, but \aL\'^/^'^^~^^^ is compa- 
rable to i?, we can no longer use series expansions to com- 
pute the binding energies. Nevertheless, we can ignore 
the effects of the non-resonant channels (with I ^ L) in 
Eq. pT|) to a good approximation. Keeping the first two 
terms in the effective range expansion for fc^-^+^ cot 5l (fc), 
and noting that k = in, we get 



1 



'-^K^ + (-1)^+1^2^+1 [1 ^ oF^^iKR)] ^ 0. (31) 



We use the above equation together with Eq. ((4]) to cal- 
culate the binding energies versus R for three fixed values 
of ol (large and positive, infinity, or large and negative) 
for each value of L. Our results are plotted in Fig. [1] 
For the s-wave resonance (L = 0), we have neglected the 
effective range ro; this is appropriate for broad s-wave 
Feshbach resonances |23l . |2J| . For p-wave and higher par- 
tial wave resonances {L > 1 ) , however, the effective range 
rL is important |24l . l33l . |3J, |37[ and we have kept it in 
our calculations. 

In Fig.[T]the power laws of Eq. (pS]) serve as asymptotes 
of the binding energies in the region R <^ laLJ^/^^^+i^. 

If tti is large and negative or equal to infinity, we have 
a single shallow bound state for each orbital magnetic 
quantum number m ranging from —L to L. It is the 
"bonding state" discussed in Sec. |TT1 and its parity is 
(T = (—1)™. If ttL is large and negative, this bound state 
disappears when the distance between the two scattering 
centers exceeds a critical value. 



Rr^ 



(XLr. 



V 2 



\aL 



1/(2L+1) 



(32) 



For the s-wave resonance, we recover the known result 
Re = |ao| |47| . For the p-wave or any higher partial wave 



Ei = 



2Mal ' 



i = 0, 

i> 1. 



(34) 



When R approaches Re from above, the energy of the 
anti-bonding state vanishes linearly: 



E 



2L+1) 
Mfl,|ai,ri,| 



L = 0, 



(2L+l)h (j^_j^^)^ 2.>1. 



(35) 



VI. A SIMPLE FORMULA FOR THE BINDING 

ENERGIES NEAR p-WAVE OR HIGHER 

PARTIAL WAVE RESONANCES 

Near a p-wave or higher partial wave resonance, if 
R, |aL|i/(2i+i) > |rL|i/(i-2^), the shallow bound states 
satisfy kR <C 1. Thus 



{-if+^n-'^+^ll + aFlUKR)] 



XLr. 



'-<T,R- 



0{R-^^+^K^), L>1, (36) 



where Uz = cr/(— 1)™ is the z-parity of the wave func- 
tion; Uz equals 1 for bonding states and —1 for anti- 
bonding states [see Eq. p^ ]. Further using the inequal- 
ity R^^^^^K^ <^ I^^lI^^ and Eq. ([5]), we can approximate 
Eq. dSIl) as 



^ _ r^ ^ Xl™^^^-2L-i ^ ifi>l. (37) 
ql n'^ 2 

From the above equation we obtain a simple formula for 
the binding energy at L > 1: 



E = -sign(aL) - <^z i j^ 



2L 



P, 



(38) 
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FIG. 1: (color online) The binding energy S of a single atom interacting resonantly with two scattering centers separated by 
a large distance R, computed using Eq. Q and Eq. (|3ip . Solid lines are for scattering volumes a_L = ±cxd. Dashed lines: a^ 
is large and negative. Dot-dashed lines: ol is large and positive. The black, red, blue, and brown curves represent orbital 
magnetic quantum numbers tti = 0, m = ±1, ra — ±2, and m — ±3, respectively. The horizontal dotted line indicates the 
binding energy Ei due to a single scattering center [see Eq. p4|l ]. We have chosen \aL\^ /\rL\^ = 50 or cxd in the 

plots for L > 1. 
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FIG. 2: (color online) The dimensionless energy E = E/\Ei\, where Ei is defined in Eq. (|34|l . versus the parameter P = 
(2L + 1)!!(2L — l)!!|a_L 1/7?^^^^ near p, d, or /-wave resonance. Curves were computed by solving Eq. (JSTJ numerically, but they 
have no visible differences from the predictions of the simple formula in Eq. (|38p . Curves with vertical intercept —1 are for 
aL > and those with (extrapolated) vertical intercept -|-1 are for a_t < 0. Curves with negative slopes are for bonding states 
and those with positive slopes are for anti-bonding states. Different colors indicate different values of ?n as in Fig. [T] in each 
graph, the curves with steeper slopes have smaller values of \m\. 



where 

P_ ^ ^_ (2L+l)!!(2L-l)!!|a,.| 
^=^' ^= i?5^T^ ■ ^^^^ 

i^i is the binding energy due to a single scattering center 
with scattering volume \aL\, and is given in Eq. (p4|) . We 
will refer to the parameter P as the "proximity parame- 
ter" . For a given separation between the two scattering 
centers, if we tunc the interaction closer to resonance, the 
proximity parameter becomes larger. For a given scatter- 
ing volume, if we place the two scattering centers closer, 
the proximity parameter gets larger. 

Equation (j38p illustrates the concept of universality in 
p-wave and higher partial wave Feshbach resonances. By 
normalizing the binding energy E by the energy due to a 
single scattering center with a positive scattering volume, 
and normalizing the distance by the scattering volume, 
we get the same set of universal relations between the di- 
mensionlcss parameters E and P regardless of the atomic 
species or the Feshbach resonance concerned. When we 
plot E against P, we shall find a set of straight lines 
with intercept — 1 for a^ < or +1 for ol > 0; the ab- 
solute values of their slopes arc given by the binomial 
2L \ _ f 2L 
L — m J \ L + m 
TO = having the steepest slopes for a given L. For ex- 
ample, for the d-wave resonance, the absolute values of 
these slopes have ratio 1:4:6:4:1. 

In Fig. [5] we re-plot the energies obtained numerically 
from Eq. (pij) . For the parameters we used [see the cap- 
tion of Fig. [1], we find that the resultant curves have no 
visible differences from the predictions of Eq. 



coefficients 



with the lines for 



VII. EFFECTS OF THE VAN DER WAALS 
POTENTIAL 

In the previous sections we assumed that the inter- 
action potential between the atom and each scattering 
center strictly vanishes beyond a certain radius r^. If 
each scattering center is a pinned neutral atom and in- 
teracts with the free atom by a long range Van der Waals 
potential 



Vir) 



2Mr6' 



r '3> re 



(40) 



Eq. ([T]) as well as the more general universal formula 
Eq. ([55]) remain valid near p-wave and d-wave resonances 
(L = 1,2), but break down for L > 3. The validity of 
Eqs. ([T]) ([55]) depends on the validity of the effective range 
expansion, Eq. (J13p . For neutral atoms with a long range 
Van der Waals potential, the effective range expansion 
remains valid up to the first two terms on the right hand 
side of Eq. p^ in the s-wave channel |48|. |49| . Near or 
at a p-wave resonance, Gao's work [4^, [5^ implies that 
the p-wave scattering length and effective range can also 
be defined [33, EE [S^ ■ Ref. [5l| discussed the conditions 
under which the effective range expansion remains valid. 



Near a d-wave resonance, Gao's work (50|, |52| also im- 
plies that the first two terms of the effective range ex- 
pansion for the d-wave phase shift are applicable within 
a certain window of collision energies. We have confirmed 
this by doing some simple perturbative expansions of the 
wave function (see Appendix Elfor details). In particular, 
we got 



tan(52(fc) = — — fc^ 1= 

^ ^ 315 900V2r(5/4)2 [j 



fcS 



ck^ 



o(fc2)]' 

(41) 

where k = k/Sg, and the parameters j and c are de- 
termined by the short-range physics. Equation Mil) is 
consistent with Gao's previous results [50, 112, 113 • F^'" 
safety we have verified Eq. (|4ip numerically by solving 
the Schrodinger equation with the —Ce/r^ potential in 
three different cases: j is zero, j is small and negative, or 
j is small and positive. According to Eq. (j^Tj) . near a d- 
wave resonance (so that j is small) and at |j| ^ |A:| ^ -, 

the term ofcfc^ is negligible, and we reproduce two terms 
of the usual d-wave effective range expansion [Eq. ([1 
if we set 



02 



TT^i 



900%/2r(5/4)2j ' 



r2 



1800V2r(5/4)2c 



(42) 



Thus, the first two terms in the effective range expan- 
sion [Eq. (J13p ] remain applicable near s,p,d- wave Fes- 
hbach resonances. Using these formulas for the phase 
shifts and the approach developed in Sec. lU we can de- 
rive Eq. ([1]) as well as the more general universal formula 
Eq. (|38|) in the presence of the Van der Waals potential 
for p-wave or d-wave resonances. 



VIII. SUMMARY AND DISCUSSIONS 

In summary, we have investigated the shallow bound 
states of a single atom interacting strongly with two scat- 
tering centers separated by a large distance R. 

At s, p, d, or /-wave resonances, we have obtained 
systematic large-i? expansions for the wave numbers 
of the bound states. Effects of non-resonant partial 
wave channels and the shape parameters in the effec- 
tive range expansion enter as correction terms in these 
expansions. For a p-wave or higher partial wave res- 
onance (i > 1) the leading term for the binding en- 
ergy behaves like — l/i?2^+i which, if combined with the 
Born-Oppenheimer approximation for two heavy and one 
light atoms, does not support three-body Efimov effect. 
Refs. [3J, Is^l contain more general discussions of the 
question of Efimov effect near higher partial wave res- 
onances, without using the BO approximation. 

Slightly off resonance we have found the 0(l/ai) cor- 
rections to the wave numbers for the bound states. 

For large scattering volumes a^ and large distances 
R we have computed the energies of the shallow bound 



states approximately, by ignoring the non-resonant par- 
tial wave scatterings and the shape parameters in the 
effective range expansion. 

Finally, for p-wave and higher partial wave resonances, 
we have found a simple formula for the binding energies, 
Eq. (|351), assuming that R and |aL|^/(^-^+^) are both large 
but may be comparable to each other. Experimental con- 
firmation of this formula will demonstrate the notion of 
universality in higher partial wave resonances. 

Although we mainly considered short-range interac- 
tions between the scattering centers and the atom, we 
have found that a long range Van der Waals potential 
will not undermine our central results, such as Eqs. ([1]) 
([55)1 . near p and d-wave resonances {L = 1,2). 

It is worthwhile to point out a technical limitation of 
our present work. We have not taken into account the 
spatial anisotropy of the interaction near a p-wave or 
higher partial wave magnetically-tuned Feshbach reso- 
nance. However, if the magnetic field is parallel to the 
line connecting the two scattering centers, our results 
remain applicable; for the shallow bound states with or- 
bital magnetic quantum number m, one will replace a;, 
n, etc in our formulas by aim and rim etc. Moreover, 
if the higher partial wave resonance between the atom 
and each scatterin g ce nter is physically a confinement- 
induced resonance [43J, |4J| due to an underlying s-wave 
Feshbach resonance, one can have a nearly isotropic effec- 
tive interaction - provided that a spherically symmetric 
trap is used to confine each pinned atom which, in turn, 
serves as a scattering center. 

To realize the scenario of one atom interacting with two 
fixed centers, it is not necessary to use two atomic species 
with a large ratio of bare masses. If two atoms of species 
A are pinned at two lattice sites of an optical lattice, 
but an atom of species B is not attracted or repelled 
by the lattice, we can treat atoms A as infinitely heavy 
scattering centers for atom B. 
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Note Added 

Near completion of this work an analogous study fo- 
cused on p-wave resonance appeared online |54| . 



Appendix A: d-w^ave phase shift in the presence of a 
Van der Waals potential 

In this appendix we derive a formula for the d-wave 
scattering phase shift 62 for a potential that behaves as 



Eq. (go]). 

At r ^ Tg the d-wave radial Schodinger equation is 

-""M + (^-7l)"W = fcV^): (Ai) 

where u(r) = rip{r). In the following we set the Van der 
Waals length /36 = 1 for convenience. 

Assuming that k <^ 1, we solve Eq. (|A1[) in two regions: 

• Region I: r ^ 1/fc; 

• Region II: r ':^ 1. 

Then, we require the two solutions to be the same in the 
intermediate region, 1 <C r ^ 1/fc. From this condition 
we obtain an approximate formula for the d-wave phase 
shift (52. 

In the region r <C 1/fc, we can treat the term k'^u{r) in 
Eq. (jAip as a perturbation, and expand u{r) in powers 
of fc^: 

n{r) = fir) + j~f{r) + fc-2[g(r) + fgir) + c~f{r)] + 0{k^) 

(A2) 
under a suitable normalization. Here 



fir)^V-rJ,„{-^) 



and 



/>) = V^J- 



5/4 



2r2 



(A3) 



(A4) 



are two independent solutions of the zero-energy radial 
Schodinger equation 



Hf{r) - 0, 



(A5) 



where H = —377 + p- — ;;i, and J„(a;) is the Bessel 
function. The functions .g(r), g{r) are the solutions to 
the following equations 



Hgir) - /(r), 
H~g{r) - Kr). 



(A6) 
(A7) 



subject to the conditions Eq. (JAIOI) and Eq. (|A11[) below. 
At large r, the functions /(r), /(r), g{r), g{r) have the 
following asymptotic expansions: 



/('•) 



5V2r(5/4) 7-2 



+ 0(r-6). 



(A8) 






30V2r(5/4) 



Kr) 



2V2 _^, 



7r(-l/4) 21\/2r(-l/4) 



(AlO) 

r + Oir"^). 

(All) 



10 



The constants j and c are determined by the short-range 
physics. At resonance, j = and the zero-energy solution 
for u{r) decays Hke r~^ at large r, consistent with the 
picture of a bound state at threshold. Slightly away from 
d-wave resonance j is close to zero and may take either 
sign. 

In the region r 3> 1, we can treat the term ^u(r) in 
Eq. ()A1|) as a perturbation, and obtain another expansion 
for u{r): 



u{r)=A{k)\ Y.k^'^Vnikr) 



n=0 



[tanS2ik)] V/c^"W„(fcr) I, (A12) 



n=0 



where A{k) is a normalization factor depending on k only, 

Vb(0-02(0, Wo{0=^n2iO, (A13) 

J2(0 E^nd 712 (^) are the spherical Besscl functions, and 
Vn{£,) and Wn{£,) are defined by the equations 



vJHO 



T-l)^"(^)-4^"-i(^)' (A14) 



e 



WL 



iO+{i-^)Wn{0 = ^Wn-l{0, (A15) 



^^ 



subject to the condition that Ki(0 ^^'^ W/n(C) a-re oscil- 
lating decaying functions at ^ — >■ (X) for n > 1. 
When ^ — > 0+ we have 



vo{o = ^+oie) 



WoiO 

^i(C) = 

w^i(e) 



_3__i_e 



0(C^ 



1 



12^& 



O 



1 



(A16) 
(A17) 
(A18) 
(A19) 



In the intermediate region, 1 ^ r ^ 1/fc, we obtain 
from Eq. ((X2|) 



w(r)cx [j + c/s2 + 0(fc4)]r3 



20\/2r(5/4)2 r2 
(other powers of r) (A20) 



and obtain from Eq. (JA12P 



u(r) oc — 4-o(fc*^) + 0(fc^tan52) 
L15 



+ { " Ss + "^^''^ + [^ + °^^^1 ^^""^'V 



(other powers of r) . 



(A21) 



Note that the "other powers of r" include terms of the 
order r"(lnr)'" for integers n and m, with n ^ —2,3. 
Comparing the above two formulas for u{r) in the inter- 
mediate region, we get the expression Eq. (|^ for the 
phase shift 62- 
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